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Abstract 



A collision of two-component Bose-Einstein condensates in the presence 
of Raman coupling is proposed and studied by numerical simulations. Raman 
transitions are found to be able to reduce collision-produced irregular excita- 
tions by forming a time-averaged attractive optical potential. Raman transi- 
tions also support a kind of dark soliton pairs in two-component Bose-Einstein 
condensates. Soliton pairs and their remnant single solitons are shown to be 
controllable by adjusting the initial relative phase between the two colliding 

condensates or the two-photon detuning of Raman transitions. 
PACS Numbers: 03.75.Fi, 32.80.Pj 
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I. INTRODUCTION 



The realization of high-density Bose-Einstein condensates invokes the study of non- 
hnear dynamics of matter waves and inventions of various optical techniques to manipulate 
Bose-Einstein condensates as well as constructions of various matter-wave interferometers 
for precision measurements. Usually, people would like to reduce irregular excitations and 
preserve the coherent property of high-density Bose-Einstein condensates as well as possible 
in many manipulating processes. However, due to the interactions between atoms, the strong 
repulsive interactions between Bose-Einstein condensates often lead to irregular excitations, 
especially in colliding processes. Inventing a convenient optical technique, which can reduce 
irregular excitations as well as exert effective controlling on the Bose-Einstein condensates, 
is very important. 

Raman coupling is a very successful experimental technique in manipulating ultra-cold 
atomic gas between different internal and external states because it has the merit of avoid- 
ing spontaneous emission loss of atoms in transition processes. For example, recently sev- 
eral experiments have already successfully employed Raman coupling in manipulating Bose- 
Einstein condensates 0. 

Here by proposing a Raman coupling scheme of atomic Bose-Einstein condensates, we 
show that Raman coupling can also be used to reduce irregular excitations produced in 
a colliding process of two Bose-Einstein condensates. Instead of irregular excitations, the 
Raman coupling produces a new type of coherent textures, namely dark soliton pairs in 
the colliding process. Additionally, we also find that these dark soliton pairs and their 
remnant dark solitons are very sensitive to the initial relative phase between Bose-Einstein 
condensates and the two-photon detuning of Raman transitions. 

The paper is organized as follows. In Sec. II, we introduce the proposed colliding process 
of two Bose-Einstein condensates in the presence of Raman coupling, and the two coupled 
Gross-Pitaevskii equations for the description of this colliding process. In Sec. HI, first we 
do some simplifications, then we show numerical simulations of the colliding process, analyze 
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the function of Raman coupling in the colhding process, and discuss the physical origin and 
some special properties of dark soliton pairs. In Sec. IV, we draw conclusions as well as 
point out some potential applications of the proposed colliding process. 

II. THEORETICAL MODEL 

First, we introduce the proposal in detail by using a schematic diagram in Fig. |[ The 
black frame, in Fig. [^(a), represents a transverse confining potential Vq{x, y) for two atomic 
Bose- Einstein condensates 1 and 2 in the abstract. This potential can be formed, for exam- 
ple, by a blue-detuned hollow laser beam propagating along z-axis, which has already been 
demonstrated to be very effective in the recent experiment Additionally, there are two 
red-detuned Gaussian laser beams 1 and 2 propagating along x-axis. The intensities of two 
laser beams are uniform in x and y directions, but in Gaussian shape in z direction, as indi- 
cated by the gray levels. Their frequencies are uji and uj2 respectively, as shown in Fig. [^(b). 
As to the condensates, we assume that atoms have three internal states, \gi)i \g2) and |e). 
Thus the two laser beams couple them together and form two Rabi transitions respectively, 
as shown in Fig. |l](b). Because of the negative detuning (A < 0) and nonuniform intensities 
of the two laser beams, they can form two z-axial confining potentials for condensates 1 and 
2 respectively, which are assumed to be initially distributed in states \gi) and \g2) respec- 
tively. Additionally, we assume that the laser beams 1 and 2 are initially separated by a 
distance and moving at velocities vi and V2 respectively in opposite directions. Thus the 
two condensates are trapped by the two laser beams and moving at velocities same as those 
of laser beams in opposite directions. While the two laser beams overlap, the condensates 
in the overlap region are transferred between states \gi) and \g2) due to the occurrence of 
Raman transitions. As the two condensates go very close, they collide and interfere with 
each other. 

Next, we derive the equation for the description of this proposal. As shown in Fig. |l|(b), 
because atoms have three internal states, first we generally consider a condensate consist- 



3 



ing of three components. To avoid spontaneous emission in Raman-transition process, we 
assume that the intermediate detuning A of laser beams 1 and 2 is much larger than other 
characteristic frequencies, such as the natural line width of the atomic transitions and the 
two- photon detuning 5. As a result, the evolution of the |e)-state component adiabatically 
follows the \gi) and |5f2)-state components, i.e., *e = -(^le"*'^^**^! Q2e~''^2**g2)/(2A), 
where ^e^^gi and are the macroscopic wave functions of the three-component Bose- 
Einstein condensate in states |e), \gi) and \g2) respectively. Then, at zero temperature and 
under rotating-wave approximation, the light-coupled \gi) and |(?2)-state components can be 
described by the following coupled Gross-Pitaevskii equations: 

(1) 

ifi-gf = + K*92 + Voix, y)^g2 + Vg2{z, t)^g2 + Uq]^ g2\^^ g2 (2) 

+C/o|*pi|'*s2 + hR*{z,t)e'^'^'-'''^'^gi 

where m is the atomic mass. Uq describes the atomic interaction strength, which is related 
to the s-wave scattering length Ugc by Uq — A.'KU^asd'^- Here we assume that age > and it 
is same for all collisions between atoms in same or different internal states. Vo(a;, y) denotes 
the transverse confining potential, which is common for both \gi) and |gi2)-state components. 
Vgx{z,t) = ^|Qi(r, i)p/(4A) and Vg2{z,t) = ^1^2(1", ^)P/(4A) are two localized optical po- 
tentials produced by the laser beams 1 and 2 respectively. R{z,t) = f2*(r, t)f22(r, t)/(4A) 
is the effective Rabi frequency of two-photon transitions. Because the laser beams 1 and 2 
are in Gaussian shape and moving along z-axis, their corresponding Rabi frequencies can 
be written as fii(r,t) = l]ioe"(^~^°i~''i*)'/"?e-^'^i"" and 1^2(r,t) = l]2oe~^^"^°^~''2*)'/"2e-^»^2--, 
where Qio and Q20 are maximum magnitudes of the two Rabi frequencies, Zqi and Z02 
are initial center positions, ai and 02 are half widths, and ki and \s.2 are wave vectors 
of the two laser beams respectively. Additionally, we assume that the wave vectors ki 
and k2 are approximately equal, then the two-photon Rabi frequency can be written as 



R{z,t) = r2iofi2oe~(^"^"^"''^*^'/''i"(^"^°'"'''*^'/"2/(4A). It is evident that when and only when 
the two potentials, Vgi{z,t) and Vg2{z,t), overlap with each other, the effective two-photon 
Rabi frequency becomes finite. All together, we reduce the description of the Bose-Einstein 
condensate from three components to two components now. 



III. NUMERICAL SIMULATION 

For simplicity and to give more prominence to the analysis of the affection of Raman 
transitions on the collision, we need to do some simplifications. We assume the transverse 
confining potential Vo{x,y) is so tight and the z-axial confining potentials, Vgi{z,t) and 
Vg2{z,t), are so loose: that in the description of the transverse modes of the Bose-Einstein 
condensate, we can neglect the affection of Vgi{z,t) and Vg2{z,t) as well as interactions 
between atoms; that in the description of the longitudinal modes, we can find that the 
longitudinal dimension of the condensate is much larger than the healing length. We can 
also find that the transverse dimensions of the condensate are very small, and the time 
scale for adjustment of the transverse profile of the condensate to the equilibrium form 
appropriate for the instantaneous number of atoms per unit length is small compared with 
the time for an excited pulse to pass a given point. This is a low-density approximation, 
in fact, similar to that used in Ref. 0. Thus by letting \I'gi(r, t) = fgi{z,t)ggi(x,y)e~'^^^^^^ 
and '^g2{r,t) = fg2{z,t)gg2{x,y)e~^^^*^^ and through a deduction similar to that in Ref. 0, 
we can get simplified one-dimensional Gross-Pitaevskii equations, 

= + + u'olUl'U + u'^lUl'U + hR{z,t)e-^^^--^^^'l'^fg2 (3) 

= + + + K\fg2?fg2 + U:\fg,\'fg2 (4) 

where Qgi and gg2 are the normalized transverse modes of the two components, /ii and /i2 
are the chemical potentials of Qgi and gg2 respectively, fgi and fg2 are the longitudinal wave 
functions, U'^ = Uq f \ggi{x, y)\^dxdy and f/g = f/o / \ggi{x, y)gg2ix, y)\^dxdy. 



As an initial condition for the colliding process, we assume that the two-component 
Bose-Einstein condensate is initially separated into two condensates, i.e., '$gi{r,t)\t=o and 
g2{r,t)\t=o, which are trapped in the ground states of the two moving potentials, Vgi{z,t) 
and Vg2{z,t), respectively. For simplicity, we assume that the two potentials have same 
shapes, depths and speeds (but in opposite directions), thus the transverse and longitudinal 
modes of the two initial ground-state Bose-Einstein condensates are in same shape and the 
condensates has same numbers of atoms and same chemical potentials. 

Then we use split operator method to solve the above time dependent Gross-Pitaevskii 
Eqs.|] and ??. First, we show a typical colliding process in the presence of Raman coupling 
in Figs. 0(a), (b) and Fig. |^. For comparison, in Figs, ^(c) and (d), we also show another 
different collision in which two condensates same as the above are trapped by two potentials 
similar to Vgi{z,t) and Vg2{z,t), but there are no Raman transitions between them. The 
elimination of Raman transitions is possible so long as the two laser beams couple the two 
ground states \gi) and \g2) with different excited states respectively instead of one common 
excited state. The parameter values used in the numerical simulation are as following: 
fi?o/(4A) = nlJ{AA) = r]io^^2o/(4A) = -6000/r, v, = -V2 = 9a,/T, a2 = a,, U', = U'^ = 
2000/i/r, zqi = —zq2 = — Oi, and the chemical potentials of initial ground-state Bose-Einstein 
condensates are both /iq = —3.99 x lO^h/r, where ai and r = 2ma{/h are considered as the 
unit length and the unit time respectively. 

In the collision without Raman transitions, the two condensates obstruct their passing 
through each other because of repulsive interactions between atoms, as shown in Figs. H(c) 
and (d). After the collision, we can see that small parts of condensates radiate out of the 
axial optical potentials, the remnant Bose-Einstein condensates oscillate violently in the 
potentials. The violent oscillations appear to be very irregular. In contrast. Figs. 0(a), (b) 
and Fig. |^ show that while Raman transitions exist, the radiations become very weak, the 
violent oscillations of remnant Bose-Einstein condensates are much reduced. This indicates 
that the repulsive interaction between two Bose-Einstein condensates is greatly counteracted 
by an attractive optical potential produced by Raman transitions. We can understand 
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the formation of this potential in this way. Because Raman transitions can transfer each 
condensate between two components, fgi and fg2 circularly, each condensate interacts with 
the two axial optical potentials Vgi and Vg2 alternatively. As a result, the light-shifted energy 
of each condensate, produced by the two laser beams, becomes a temporal combination of 
Vgi and Vg2, and each condensate feels a vibrating optical potential locating at the collision 
center z = 0. According to the overlap form of Vgi and Vg2, the vibrating potential, in 
average, is attractive. Thus this time-averaged attractive potential counteracts the repulsive 
interaction between two condensates, so it becomes much easier for them to pass through 
each other. We therefore observe the oscillation amplitudes of condensates, produced by the 
repulsive interaction, become much small. But between t = 0.05r and 0.15r, we can still 
observe that the volume of each condensate oscillates shghtly in the colliding process due to 
the vibration of the attractive optical potential. 

In the colliding process, we also observe that dark solitons are generated, as shown in 
Fig. 1^. Formations of dark solitons are due to the interference of colliding condensates. 
Although the two condensates are initially distributed in two different internal states re- 
spectively, Raman transitions can transfer parts of them from one of the internal states to 
the other in the colliding process. Consequently, when the transferred part of one conden- 
sate overlaps with the part of the other condensate in the same internal state, interference 
fringes are produced. Due to nonlinear effect of condensates, interference fringes evolve into 
dark solitons. It is evident that this process is similar to that in Ref. 0. 

However, solitons produced in this colliding process in the presence of Raman coupling 
are quite different from those single solitons described in Ref. One of important 

differences is that these solitons always appear in pairs. As shown in Figs. |^(a)-(d), in the 
colliding process, every soliton in one component of the condensates always has a corre- 
sponding soliton in the other component at the same spatial position. Another important 
difference is that the density difference between two components does not lead to a speed 
difference between two solitons in a pair. As shown in Figs.|^(a) and (b), sometimes the 
two components have quite different local densities at the location of a soliton pair, how- 
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ever, solitons in the pair always move synchronously in the condensates. The reason for 
synchronous motion of two solitons in a pair is that Raman transitions associate the two 
components tightly and consequently locks their phases together. As we can see in Fig. ^, 
the phases of two components are almost completely same in the region of strong Raman 
transitions that are indicated by magnitudes of the normalized two-photon Rabi frequency 
|4i?A/ (i7iofi2o)|- Dark solitons are phase kinks in essence, so there is no doubt that solitons, 
with similar spatial phases, in a pair move synchronously. Additionally, the cross phase 
modulation between dark solitons in a pair also makes a repulsive effect between them, 
and tends to separate them. It is evident that this repulsive effect is also counteracted by 
Raman transitions, which produce an equivalently attractive effect between solitons in the 
pair by locking their phases together. Even after a collision of two soliton pairs A and B, 
solitons in each pair can still preserve their waveforms and trajectories well, as shown in 
Figs, ^(b)-(c) and Figs. H(a)-(b) respectively. This illustrates that collisions between soliton 
pairs are elastic collisions, which are usually considered as an important common character- 
istic of all kinds of solitons. The synchronous motion of solitons in the pair also indicates 
the propagation speed of the soliton pair is determined by the total local density of the two 
components instead of each one component. We can find this property from Figs, ^a)-(b). 
The propagation speed of the soliton pair, in the region 2; < 0, at time t = 0.078r is larger 
than that at time t = 0.047r, because the total density of the two components is increased 
in the overlap region,. Some similarities between a soliton pair and a single soliton can 
also be seen. For example, the propagation speed, the width and the density contrast of a 
soliton pair also depend on the value of its phase kink that is also varied by the local density 
gradient. We can find these properties by a careful observation of Figs. 0(a)-(b), Fig. ^ and 
Fig. ^ Because the density of the condensates varies in parabola-like form in space, we can 
find its phase kink flips in the boundary region of the condensates, simultaneously its speed 
flips too, and as a result, the soliton pair oscillates inside the condensates. 

The two components of the condensates begin to separate very apparently after t = 0.15r, 
because of the separation of the two axial optical potentials, as shown in Figs, ^a)-(b). 
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Simultaneously, Raman transitions become weaker and weaker because of the decrement 
of the overlap of the two laser beams, as shown in Figs, ^(d)-(f). As a result, soliton 
pairs become unstable and are to be separated in this process. Critically depending on the 
local evolution of the two components, some solitons disappear and some of them remain 
and propagate into the separated one-component condensates, as shown in Figs. 0(a)- (b) 
and Figs. |(e)-(f). Because of the disappearance of cross phase modulation in each one- 
component condensate, the oscillation behaviors of remnant single solitons are different 
from those of soliton pairs. Especially, the oscillation periods of remnant single solitons 
appear to be quite different from those of previous soliton pairs. 

The generation of dark soliton pairs is controllable. Usually, the relative phase between 
two condensates in different internal states is meaningless because of the orthogonality of the 
internal states. However, here Raman transitions can form a circulation of atoms between 
the two internal states, the relative phase between the two condensates are related and 
therefore becomes very important. Because the formation of dark soliton pairs is due to the 
interference of colliding condensates, we can also control the generation of soliton pairs as 
well as remnant single solitons by adjusting the initial relative phase between two colliding 
condensates or the two-photon detuning of Raman transitions, as shown in Fig. ^. 

In the above, we have discussed the collision of two Bose-Einstein condensates under 
the one-dimensional approximation, then will a three dimensional collision be completely 
different from the one- dimensional collision? Some difference must exist. For example, in 
three-dimensional process, dark soliton pairs might be distorted by transverse perturbations 
and evolve into vortices because similar behaviors of single dark solitons have been found 
P). However, in the above analysis, we have already found that the function of Raman 
coupling in the colliding process is, in fact, to lock the phases of the two components of the 
condensates by forming the strong atom circulation between the two internal states, and as 
a result, solitons appear in pairs. As for vortices, they are phase singular points, so they 
are quit similar to dark solitons in essence. We think that Raman coupling can also lock 
the phases of two components of the condensates in a three-dimensional case, and vortex 
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pairs might therefore be formed. Certainly, more affirmative answer should be gotten in a 
numerical simulation of a three-dimensional collision, which will be done in the future work. 

IV. CONCLUSION 

In conclusion, we have proposed a new Raman coupling scheme of Bose-Einstein con- 
densates, and analyzed the influence of Raman transitions on the collision and excitation 
of two-component Bose-Einstein condensates. We have found that Raman transitions can 
reduce collision-produced irregular excitations by forming a time-averaged attractive poten- 
tial. Raman transitions also support a new kind of dark solitons, i.e., dark soliton pairs in 
two-component Bose-Einstein condensates, by locking the phases of the two components. 
We have also shown the control of soliton pairs and their remnant single solitons by ad- 
justing the initial relative phase between the two colliding condensates and the two-photon 
detuning of Raman transitions. The dynamical sensitivity of the dark soliton pairs and 
remnant solitons to these parameters indicates that the proposed collision can be used in 
the measurement of relative phase between different high-density Bose-Einstein condensates. 
Currently, people are considering constructions of various interferometers with matter waves 
for high precision measurements. We think that this proposal might provide a useful method 
to this area. 
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FIGURES 

FIG. 1. (a) Collision and Raman coupling scheme of two atomic Bose-Einstein condensates, 
(b) Atomic internal energy levels of Bose-Einstein condensates and light coupling, j^i), \g2) denote 
two ground states and |e) denotes an excited state, energy eigen values of which are 0, hujg and fiuje 
respectively. The laser beams 1 and 2, with frequencies uji and L02 respectively, couple these states 
together, forming so-called "A" configuration Raman transitions, ili and ^2 are corresponding 
single-photon Rabi frequencies. The intermediate detuning and two-photon detuning of Raman 
transitions are A = — and 5 = loi — 0J2 — i^g respectively. Direct transitions between \gi) and 
\g2) are electric-dipole forbidden. 

FIG. 2. Contour plots of time evolution of the condensate wave functions \ fgi\ and \ fg2\ in two 
colliding processes, (a) is for \fgi\ and (b) is for \fg2\ in the collision with Raman coupling, (c) is 
for \fgi\ and (d) is for \fg2\ in the collision without Raman coupling. In these collisions, the initial 
phases of fgi and fg2 are 7r/4 and respectively, and the detuning 5 = 0. The dash-dot lines denote 
the boundaries of the corresponding axial optical potential, Vgi{z,t) or Vg2{z,t). These boundaries 
correspond to the boundaries of Gaussian laser beams at the 1/e of maximum intensities. Between 
t = 0.05r and 0.15r, there are four zigzag density dip traces in (a) and (b). These four density dip 
traces indicate the oscillations of four dark solitons inside the condensates. 

FIG. 3. Generation of dark soliton pairs in the collision with Raman coupling and remnant 
single solitons. The thin solid lines denote \fgi/fo\, the thick solid lines denote \fg2/fo\, and the 
dashed lines denote \4:RA/{^liQfl2o)\- fo is the maximum of |/gi|t=o and \ fg2\t=o- There are four 
dark solitons in each figure of (a), (b) and (c). The arrow A marks one soliton pair, and the arrow 
B marks another soliton pair respectively. The values of t denote the time of subfigures sampled 
from the colliding process shown in Fig.||(a) and (b). 
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FIG. 4. Phase evolution of the condensate wave functions in the colhsion with Raman coupling. 
The solid lines in (a)-(c) denote phase of fgi, the solid lines in (d)-(f) denote phase of fg2/ fgi, and 
the dashed lines denote \4RA/{0,iqQ2o)\- The arrow A marks the position of one soliton pair, and 
the arrow B marks the position of the other dark soliton pair respectively. The values of t denote 
the time of subfigures sampled from the colliding process shown in Fig.|2|(a) and (b). 

FIG. 5. Contour plots of time evolution of the condensate wave functions \fgi\ and \fg2\ in two 
colliding processes, (a) is for \fgi\ and (b) is for \fg2\ in one collision. In this collision, the initial 
phases of fgi and fg2 are 7r/2 and respectively, and the detuning 5 = 0. (c) is for \fgi\ and (d) is 
for \fg2\ in the other collision. In this collision, the initial phases of both fgi and fg2 are Os, and the 
detuning 5 = 100/r. The dash-dot lines denote the boundaries of the corresponding axial optical 
potential, Vgi{z,t) or Vg2{z,t). These boundaries correspond to the boundaries of Gaussian laser 
beams at the 1/e of maximum intensities. 
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